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A topological phase can often be represented by a corresponding wavefunction (exact eigenstate of a model 
Hamiltonian) that has a higher underlying symmetry than necessary. When the symmetry is explicitly broken 
in the Hamiltonian, the model wavefunction fails to account for the change due to the lack of a variational 
parameter. Here we exemplify the case by an integer quantum Hall system with anisotropic interaction. We 
demonstrate the recovery of the variational parameter in a single-mode approximation, which is consistent with 
the recently proposed geometric consideration of the quantum Hall phases. 



Introduction. The integer quantum Hall (IQH) effect is now 
regarded as a simple class of topological insulators with bro- 
ken time-reversal symmetry. A system in such a phase has 
invariant properties of its ground state wavefunction under 
smooth deformations of the Hamiltonian. This, of course, 
does not mean that the ground state wavefunction itself is in- 
variant. For convenience, we often choose to study the system 
with continuous symmetries in translation and/or rotation in 
which the ground state wavefunction is relatively simple, al- 
though these symmetries are unnecessary for the study of the 
topological properties. In fact, such a choice can sometimes 
hide important physics. 

Consider an IQH system in the disk geometry with ro- 
tational symmetry. Electrons in a polarized lowest Landau 
level (LLL) fill up ring-shaped orbitals 4>m{z) = 2™e~l^l /^, 
where z = x-\-iy\s the complex coordinate on the plane. The 
trivial many-body ground state is the Slater determinant of the 
orbitals with the lowest possible m 
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which is the well known Vandermonde determinant multiplied 
by the ubiquitous Gaussian factor The wavefunction, which 
vanishes when two electrons coincide with each other, man- 
ifests the Pauli exclusion principle. The trivial many-body 
state is known to be the ground state even in the presence of 
interaction. 

This prompts an interesting question: Can effects of inter- 
action be explicitly manifested in the robust IQH state? If yes, 
how do they modify the ground state wavefunction variation- 
ally? In fact, to construct a variational wavefunction for cor- 
related systems from its noninteracting counterpart is a long- 
standing problem uJ]. Suppose we decompose a many -body 
Hamiltonian tohe H — Hq + Hi. We assume that we know 
the ground state \^q) of Hq, which describes, for the present 
discussion, the noninteracting part of the total Hamiltonian. 
Hi describes the interparticle interaction. Within the single- 
mode approximation, which replaces the effect of Hi by a sin- 
gle, mean excitation energy ujq, we can write down an ansatz 
wavefunction |^o) — e~^i/'^°|$o) for the ground state of 



H. Here, luq can be tuned as a variational parameter to min- 
imize the ground state energy. A well-known example is the 
Gutzwiller's wavefunction yj, which accounts for the elec- 
tron correlation in the on-site Hubbard model. An improved 
independent-mode approximation can lead to Jastrow's ansatz 
for a trial ground-state wavefunction 

*o (ri, ..., tn) = e^- ■^('■'-'■^■)$o (ri, ..., r^) , (2) 

where the Jastrow function /(r — r') is determined by energy 
minimization |l|]. According to this description, the Jastrow 
function for the IQH state simply vanishes despite the inter- 
particle interaction. 

This, of course, is counterintuitive. In fact, the absence 
of the Jastrow factor for the isotropic IQH state, which can 
be regarded as a special case of the Laughlin state at fill- 
ing fraction i^ = 1, has already been questioned in the con- 
text of the geometric description of fractional quantum Hall 
(FQH) states [4, 5]. Haldane [4] argued that guiding center 
metric should be regarded as a variational parameter in the 
generalized family of Laughlin trial wavefunctions. Explicit 
construction of the family of wavefunctions ySl] confirms that 
the isotropic Laughlin wavefunction is the variational ground 
state for an isotropic interparticle interaction. In contrast, the 
dipole-dipole interaction, as occurs in rotating dipolar cold 
atomic systems, leads to an anisotropic variational ground 
state 1 5]. On a torus geometry without boundary, one can cali- 
brate the wavefunction anisotropy or the guiding center metric 
by ejcploring static structure factor and pair correlation func- 
tion illl. 

The geometric description of quantum Hall states echoes 
the single-mode (independent-mode) approximation in that 
both quest for a (set of) variational parameter(s) that accounts 
for the correlation of electrons. It is, therefore, conceivable 
that for the IQH effect the two descriptions may be unified. 
It should be pointed out that the variational parameter also 
emerged from the discussion of the stress response to an ap- 
plied field by Read and Rezayi |8] for two-dimensional elec- 
trons in a uniform magnetic field. Although the authors there 
turned to reveal the illuminating relation of the Hall viscos- 
ity to the mean orbital spin per particle, their metric-deformed 
single-particle basis is sufficient to set up the stage for a van- 
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FIG. 1: (Color online) The mean total angular momentum of the 
ground state 17 = (*(^'(a, e)\L^ \'^^^\a, 6»)> (in units of h) in a 
system with TV = 6 particles for various confining potential strength 
a and tilt angle 9. The blue region to the top and right of the orange 
dashed line represents the anisotropic IQH state, while the rest is 
loosely referred to as the FQH regime, with the brown region in the 
lower left corner representing the v = 1/3 Laughlin state. The IQH 
regime is separated into the perturbative and nonperturbative IQH 
regimes, depending on whether Eq. (|9) is a good description of the 
ground state or not. Numerical analysis suggests that M — Mo =0.5 
can be used to define the boundary between the two IQH regimes. 



(3) 



ational IQH wavefunction |15| 



where 7 parametrizes the metric change or anisotropy. This 
apparently differs from Eq. (|2), but we will argue they are 
consistent, at least for the IQH state with dipole-dipole inter- 
action. 

In this Letter we take a pedagogical approach to construct 
the variational ground-state wavefunction from numerical di- 
agonalization of a microscopic system with dipole-dipole in- 
teraction. We find that in the perturbative regime the realis- 
tic ground state can be written as the isotropic IQH ground 
state multiplied by the sum of the leading terms in a geomet- 
ric series, which can be cast in an exponential form much like 
the Jastrow factor in Eq. (|2). The Jastrow factor can be un- 
derstood as arising, in the single-mode approximation, from 
an anisotropic quadrupolar interaction. Up to a center-of- 
mass contribution not generated by the dipole-dipole interac- 
tion, we show that the wavefunction agrees with the Bogoli- 
ubov transformation of the guiding center coordinates in the 
isotropic wavefunction, or Eq. ( 3}. 

Model. As proposed in Ref. |9l the fast rotating quasi-two- 
dimensional gas of polarized fermionic dipoles serves as an 
ideal arena for the present study. The Larmor theorem states 
that the fast rotation is equivalent to a high magnetic field for 
the fermions. Since the p-wave interaction for the polarized 
fermions is typically small unless in the resonance regime, the 
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FIG. 2: (Color online) Low-energy excitations for the isotropic IQH 
state with a = 0.1 and 6 = Q° . The blue thick levels can be identi- 
fied as $0 [Eq. Q], $i [Eq. Q], and $1 [Eq. ^]. Note that their 
excitation energies /S.E scale with their momenta AM. These states 
are the main components of the anisotropic IQH ground state for 



only significant interaction in the system is the dipole-dipole 
interaction. Therefore, we set up a quantum Hall system with 
anisotropic dipole-dipole interaction as described by the fol- 
lowing Hamiltonian 



H = aV 



where /^^ 

Tlim''^fmfm IS the z-componcnt of the total angular mo- 
mentum. The details of the anisotropic dipole-dipole interac- 
tion matrix elements Vi234(6') are explained in Ref. [9]. The 
isotropic confining potential strength a can be tuned by the 
rotation frequency and the tilt angle 9 by applied electric and 
magnetic field. We explore, by exact diagonalization, the A^- 
particle ground state and compute its mean total angular mo- 
mentum (in units of h) 
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M= U'^^\a,e) V *W(a,6') 
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where $*^^'(a, 6) denotes the ground state wave function for 
given parameters a and 6. 

As illustrated in Fig. [T] with isotropic dipole-dipole inter- 
action in the plane of motion, i.e. for Q = 0°, the total angular 
moment is a good quantum number and can be matched to 
that of the IQH state Mq = N(N - l)/2 for a > 0.085 
and Ml = 'SN{N - l)/2 the v = 1/3 LaughHn state for 
a < 0.018. Note that we only show a > 0.01 to avoid the un- 
wanted complication due to weak confinement and finite sys- 
tem size. As 9 increases, the rotational symmetry is explicitly 
broken and M increases from Mq to reflect the anisotropy of 
the interaction. We identify the extended blue region to the 
up and right of the curved orange dashed line in Fig. [Uto be 
the anisotropic IQH phase with M ~ Mq + 0(1), which 
has the maximum density in the interior of the oval-shaped 
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FIG. 3: (Color online) Overlap analysis of 1 - K^'^^Ha- ^)l'I'o)|^ 
|{*('^)(a,6l)|<E>i)p, and |(*(^'(q,6I)|$1)P for (a) various 6» at a 
fixed a — 0.1 and (b) various a at a fixed 6 = 70°. The overlap 
analysis suggests that "^^^'{a = 0.1, 6) is mostly composed of $o, 
<I>i, and $1. In (b) we also plot j(*''"(a, 6)\^i)f calculated from 
Eq. dlOb . which only deviates from the direct overlap calculation for 
weak confinement a < 0.03. 



droplet [9]. As labeled in Fig.lT] the IQH phase can be further 
separated into perturbative and nonperturbative regimes to be 
explained later The rest of the parameter space is loosely re- 
ferred to as the FQH regime ||9|] . 

Ground-state wavefunction. The isotropic IQH ground 
state (in a harmonic trap) is the maximum density droplet of 
electrons in the LLL. One expects that in the presence of small 
anisotropic interaction the IQH state can only be perturbed at 
the edge. Hence the anisotropic IQH state is expected to be the 
superposition of the isotropic IQH ground state and its edge 
states, whose spectrum is shown in Fig.|2]for a = 0.1. Due 
to the nature of the dipole-dipole interaction, only edge states 
with an angular momentum difference from the ground state 
by an integral multiple of 2 can be generated. Therefore, the 
leading correction must come from either one of the following 
states or both. 
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K$(«)(o.i,0)|*o)l 
l(^(«)(o.i,0)|*i)| 



0° 



20° 



40° 



60° 



80° 



1.0000 0.9994 0.9945 0.9865 0.9809 
0.0000 0.0342 0.1042 0.1628 0.1927 
0.0000 0.0008 0.0074 0.0183 0.0258 



TABLE I: Overlaps of the ground state ^^'^^ (a = 0.1,0) obtained 
by exact diagonalization of an A'^ = 6 system and the wavefunctions 



$0, '^\, and $|. 



highlight a third energy level, which corresponds to 
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n(^,-z,)e-5:.l-lV4. (8) 
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We should emphasize that the labeled edge excitations in 
Fig.|2]for realistic interaction are not exactly the model wave- 
functions, i.e., specified by Eqs. (Q and ([8]), respectively, but 
the corresponding overlaps are close to unity. For example, in 
an iV = 6 system with isotropic dipole-dipole interaction, the 
overlap of the corresponding state with $4 is 99.85%. From 
the form of the wavefunctions, as well as the linear energy- 
momentum relation, of the states <I>2 and $4, the latter is 
clearly the corresponding two-boson excitation of the former. 

First, we focus on the regime with M — Mq ^ 1, in which 
these edge-state wavefunctions are sufficient for the analy- 
sis of anisotropic IQH states. We choose a = 0.1 (the top 
boundary in Fig. [U and calculate the overlaps of the ground 
state ^(^)(a, (?) with A^ = 6 obtained by exact diagonal- 
ization and the wavefunctions $o> ^% ^nd $|. In Figure 
Ha) we plot 1- |(*(6)(a,6i)|$o)l^ \{-^^^\a,e)\^l)\^, wA 
\{^^^\a,6)\^'\)\^ as a function of 9. For convenience, we 
also list the values of overlap at selected angles 9 in Table U 
As 9 increases, the ground state evolves from $0 to a weighted 
sum of <i>o and <i>2, and eventually to a weighted sum of ^q, 
<i>2, and <i>4. At a = 0.1 the three states together exhaust 
*(^)(a, 9), e.g. up to 99.996% at 61 == 80°. 

Not surprisingly, the overlap with the higher-momentum 
<i>| is much smaller than that with $|. A question arises 
whether the former can be related to the latter in the perturb a- 
tive regime. Given the general knowledge of trial wavefunc- 
tions for correlated systems, we attemped to fit the numerical 
wavefunction by 



i<j 
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which are unnormalized but orthogonal to each other in the 
thermodynamic limit. The nomralization factor can be worked 
out using the algebra of symmetric polynomials. Numerical 
analysis reveals that the ground state has vanishingly small 
overlap with $2' but substantially large overlap with $2- I" 
fact, both $2 ^nd $0 are known to be the eigenstates of the 
projected isotropic interaction in the lowest Landau level lllOll . 
which are indicated in the low-lying energy levels for the cor- 
responding isotropic interaction in Fig.|2] In the figure we also 



of which the exponential factor can be expanded to include 
(i>o, ^2^ and $1 as the leading terms. If this is a good trial 
wavefunction, we expect 
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For A^ = 6 the right side can be evaluated to be 0.0009, 
0.0079, 0.0193, and 0.0270 for 9 = 20°, 40°, 60°, and 



80°, respectively. The calculated values are in good agree- 
ment with the tabulated ones in Table |l] This regime, with 
Eq. (|9]l as a good variational wavefunction, is dubbed as the 
perturbative IQH regime. We point out that, in the pertur- 
bative regime, the overlap with the three-boson excitation 



ubov) transformation 



n 



Ei<j(-2^» - ZjYf^o is vanishingly small. 



Figure [3tb) repeats the overlap analysis for various a 
at a fixed 9 = 70°. For comparison, we also plot 
|(*(6)(q,^61)|$2|j|2 calculated from Eq. ^. We find that 
'^^'^'{a,9) is still, to a good approximation, exhausted by 
$0, *2' and<l>|. However, \{'i!^^\a,6)\'^'l)\'^ calculated from 
Eq. ( [Tol l and by the overlap analysis deviate from each other 
for a < 0.03. The onset of the deviation is almost indepen- 
dent of 9 in the anisotropic IQH phase and coincides well with 
A/ — Mo = 0.5. Therefore, we dubbed the regime with larger 
M the nonperturbative IQH regime, as illustrated in Fig.[Tl in 
which Eq. (|9]l is insufficient in describing the true ground state 
precisely. 

The form of Eq. (|9j exemplifies the independent-mode ap- 
proximation in Eq. (|2|i in the IQH context. In particular, since 
$4 is the two-boson excitation of $|, the wavefunction Eq. (|9|l 
is a single-mode approximation for the anisotropic IQH state. 
As discussed in the introduction to the single-mode approx- 
imation, the wavefunction Eq. (|9]l suggests one can extract, 
under the single-mode approximation, a quadrupolar interac- 
tion V{zi, Zj) = {zi — Zj)^ from the anisotropic interaction. 
More formal formulations require us to write down a real in- 
teraction V{zi,Zj) + V{zi,Zj), whose antiholomorphic con- 
tribution drops out in the LLL projection. Not surprisingly, 
the real form faithfully reproduces the angular dependence of 
the anisotropic component of the dipole-dipole interaction. 

Discussions. The consideration of the broken rotational in- 
variance in the quantum Hall context appeared first in the al- 
ternative discussion on the effect of anisotropic effective mass 
tensor [11]. Very recently the joint effects of both anisotropic 
mass and interaction surfaced in the geometric description of 
quantum Hall states |4], of which the emergence of the Jas- 
trow factor can be understood as a special example. Hal- 
dane ij] pointed out that the conventional understanding of 
the Laughlin wavefunctions is not complete; in his novel geo- 
metric description the original Laughlin wavefunction is sim- 
ply a member of a family of Laughlin states, parameterized by 
a hidden (continuous) geometric degree of freedom. The fam- 
ily of the Laughlin states, with the geometric factor as a vari- 
ational parameter, potentially provides a better description of 
the FQH effect in the presence of either anisotropic effective 
mass or anisotropic interaction, which are present in real ma- 
terials. The present work also confirms that even in the IQH 
context the geometric description is necessary and we uncov- 
ered the variational parameter 7 by our numerical anatomy of 
the ground state wavefunction. 

The family of the variational wavefunction can be con- 
structed alternatively by squeezing the LL orbitals. In the 
symmetric gauge one can introduce a unimodular (or Bogli- 




(11) 



of the LL wavefunctions to generate a set of anisotropic 
single-particle basis states H. As a result, we find that 
the z^ = 1 IQH state (in the absence of LL mixing) be- 
comes Eq. (O. Note that the extra term in the exponential, 
which arises from the modification of single-particle basis, 
will be ubiquitous for all FQH wavefunctions. The unimodu- 
lar transformation approach also appeared in the consideration 
of Read and Rezayi [8] on the geometric aspects of the uni- 
form deformation of quantum Hall states in the context of Hall 
viscocity, which is the nondissipative transport coefficient that 
describes the stress in the adiabatic response to a strain. The 
authors [8J pointed out that the extra term in the exponen- 
tial corresponds to a quadrupolar harmonic perturbation in the 
plasma mapping for the Laughlin state fl2|]. 

Apparently, Eq. (O differs from Eq. (|9}. But the two forms 
can be related by noticing 
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The difference e"^^^'^'-' ' , when expanded for small 7, in- 
volves $2 and its descendents, which are not generated by the 
dipole-dipole interaction, as we discussed in the wavefunction 
analysis. Therefore, the unimodularly transformed wavefunc- 
tion for the IQH case is consistent with the single-mode ap- 
proximation of the anisotropic quadrupolar interaction. 

The explicit construction of the family of variational wave- 
functions in the FQH case was discussed by two of the authors 
and collaborators. For concreteness, the anisotropic Langhlin 
state atiy = 1/3 is 



nh 



+ 27* (a,. -a, J 
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(13) 
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where we omit the usual Gaussian factor for simplicity; we 
use d to emphasize that the partial derivatives are not oper- 
ated on the (omitted) Gaussian factor The wavefunction is 
simply the unimodular transformation of the isotropic Laugh- 
lin wavefunction projected to the LLL. For small 7 the ex- 
ponential, shared by the IQH wavefunction, generates edge 
excitations, i.e. density deformation at the edge. The Jastrow 
product describes the deformation of the correlation between 
particles, i.e. the shape of the correlation hole. As a result of 
the unimodular transformation construction, the two kinds of 
deformations are controlled by the same (complex) anisotropy 
parameter For a realistic edge, the confining potential is ex- 
pected to squeeze the shape of the correlation hole locally near 
the edge, which may offers an alternative understanding of 
FQH edge physics in the geometric description 1113 1. 

In summary, we find that the anisotropic IQH phase due 
to interparticle dipole-dipole interaction can be separated into 



the perturbative and nonperturbative regimes, based on the to- 
tal angular momentum deviation from its isotropic value. In 
the perturbative regime, we show that a variational wavefunc- 
tion can be constructed in the single-mode approximation as 
in the generic treatment for correlated systems. We further 
argue that the single-mode approximation for the anisotropic 
interaction is consistent with the geometric description for the 
anisotropic IQH phase. 
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